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ABSTRACT 

We report the strategy and results of our radial velocity follow-up campaign with the SOPHIE spectrograph (1.93-m OHP) of four 
transiting planetary candidates discovered by the Kepler space mission. We discuss the selection of the candidates KOI-428, KOI-410, 
KOI-552, and KOI-423. KOI-428 was established as a hot Jupiter transiting the largest and the most evolved star discovered so far 
and is described by Santerne et al. (201 la). KOI-410 does not present radial velocity change greater than 120 m s _I , which allows us 
to exclude at 3<x a transiting companion heavier than 3.4 Mj up . KOI-552b appears to be a transiting low-mass star with a mass ratio of 
0.15. KOI-423b is a new transiting companion in the overlapping region between massive planets and brown dwarfs. With a radius of 
1.22 ± 0.1 1 R Jup and a mass of 18.0 ± 0.92 M Jup , KOI-423b is orbiting an F7IV star with a period of 21.0874 ± 0.0002 days and an 
eccentricity of 0. 12±0.02. From the four selected Kepler candidates, at least three of them have a Jupiter-size transiting companion, 
but two of them are not in the mass domain of Jupiter-like planets. KOI-423b and KOI-522b are members of a growing population 
of known massive companions orbiting close to an F-type star. This population currently appears to be absent around G-type stars, 
possibly due to their rapid braking and the engulfment of their companions by tidal decay. 

Key words, planetary systems - brown dwarfs - binaries: eclipsing - techniques: photometry - techniques: radial velocities - tech- 
niques: spectroscopic 
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*0 . 1. Introduction KOI-428 was revealed as a transiting hot Jupiter with a radius of 

£j' T , , • . onnn v , ■ . • • 1.17+0.04 R Jup and amass of 2.2+0.4 M Jup orbiting the largest 

(N ■ Launched in March 2009, Kepler is the second space mission and ^ mQst e P yolved star (2.13 R and 1 .48 M ) discovered so 

, designed to find transiting exoplanets by high-accuracy photom- faf ^ a tranfjiti exoplanet. This planet KOI-428b, which also 

• etr y- Jhis mission has already demonstrated its strong capabil- c onds t0 the first K ler lanet established from the public 

ity of detecting transiting exoplanet candidates. Thanks to ra- d wag described b Santerne et al (201 q 

y—{ dial velocity follow-up or/and transit timing variation it has es- 

t-H . tablished and characterized up to eight new transiting planetary 

systems from the super-Earth regime - Kepler-lOb (Batalha et 2. Selection Of Kepler candidates 

. al. 201 1) - to the massive and bloated Jupiter regime - Kepler- 

X ■ 5b (Koch et al. M). From the 156,000 stars continuously ob- first set of ^6Kepler exoplanet candidates was published by 

^ • served by Kepler, 706 exoplanets candidates were first identified Bomckl et aL These candidates are faint stars (14 < 

«J ; and 306 of them were published by Borucki et al. (201 la) in m v ^ 16 ) that are not lncluded ln the ground-based observation 

February 2010. To establish the planetary nature of these candi- f°How-up conducted by the Kepler team. They were identified 

dates, to assess the fraction of false positives, and to characterize from the first 33 ' 5 " da y se 8 ment of the sclence operations (Ql) 

the true mass and density of new transiting planets, we started from Ma ? 13 to June 15 ' 2009 wlth a temporal sampling of 29.4 

a large program of radial velocity follow-up with the SOPHIE min ' 

high-resolution spectrograph mounted on the 1.93-m telescope From thls llst of candidates, we defined a selection using 

at the Observatoire de Haute-Provence Observatory (Bouchy et the following criteria. We first removed all targets fainter than 

al.[2p09 Perruchot et al. 2008). magnitude 14.7, which is close to the magnitude limit of the 

We present here the strategy and results of our first campaign S0PHIE spectrograph. One hundred fifteen candidates have a 

during third quarter 2010 on the four Kepler object of inter- Ke P ler magnitude in the range 13.9 - 14.7. We then extracted 

est KOI-428, KOI-423, KOI-552, and KOI-410. Our first target candidates with an estimated radius greater than 0.8 R Jup . Almost 

all the known transiting planets with R > 0.8 Rj up have a mass 



Send offprint requests to: F. Bouchy e-mail: bouchy@iap . fr higher than 0.4 M Jup . Lower masses are difficult if not impos- 

* Based on observations made with the 1.93-m telescope at the sible to detect with SO PHIE in this range of stellar magnitude 
Observatoire de Haute-Provence (CNRS), France (Santerne et al. 1201 lbi l. We removed the monotransit candidates 
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Table 1. Kepler candidates from Borucki et al. (201 la) selected for the SOPHIE follow-up. 



Kepler object 


Kepler input catalog 


RA 


DEC 


Kepler 


Rp" 


Period 


Terr 


SOPHIE follow-up 


of interest 








mag 


[Rjup] 


[days] 


[K] 


status 


KOI-4 10.01 


KIC-5449777 


19h28m59.53s 


+40d41'45.8" 


14.5 


1.07 


7.217 


5968 


planet or DEB 


KOI-423.01 


KIC-9478990 


19h47m50.46s 


+46d02'03.5" 


14.3 


0.94 


21.087 


5992 


brown-dwarf 


KOI-428.01 


KIC-10418224 


19h47ml5.29s 


+47d31'35.8" 


14.6 


1.04 


6.873 


6127 


planet 


KOI-552.01 


KIC-5122112 


19h49m35.58s 


+40dl3'45.1" 


14.7 


1.00 


3.055 


6018 


low-mass star 



" Estimated radius of the transiting candidates 



(with an estimated period longer than 33.5 days). This led to a 
list of 1 1 candidates. The last step was to try to privilege massive 
hot- Jupiter and brown-dwarf candidates with the aims of better 
understanding this class of massive objects, of further exploring 
the link between massive planets and brown dwarfs and of de- 
riving the main properties of their parent stars. As discussed in 
Bouchy et al. (2011), almost all transiting massive planets and 
brown dwarfs are found to orbit F-type stars as a consequence of 
rapid rotators, independent of their age (Nordstrom et al. 1 19971 ). 
We then selected the four hottest stars in our list for our first cam- 
paign with an estimated effective temperature greater than 5900 
K. The characteristics of these four targets are listed in Table [TJ 
More recently, Borucki et al. (201 lb) have published an up- 
date of the Kepler exoplanet candidates from the analysis from 
the analysis of the second quarter (Q2) data lasting 88.7 days 
from June 20 to September 16, 2009. We thus downloaded from 
the MAST databas43 these two first quarters of the publicly 
available Kepler data for our selected candidates. We then com- 
bined our radial velocity observations to the analysis of the can- 
didates' Kepler light curves. 

3. SOPHIE observations 

We performed the radial velocity (hereafter RV) follow-up of the 
four candidates with the SOPHIE spectrograph (Bouchy et al. 
2009 Perruchot et al. 2008 ) installed on the 1 .93m telescope in 
the Observatoire de Haute-Provence, France. SOPHIE is a cross- 
dispersed, high-resolution fiber-fed echelle spectrograph, stabi- 
lized in pressure and temperature and calibrated with a thorium- 
argon lamp. It is mainly dedicated to measure high-accuracy ra- 
dial velocities on solar-type stars for exoplanets and asteroseis- 
mology studies. SOPHIE is one of the key facilities in the radial 
velocity follow-up of CoRoT and SWASP-North transiting plan- 
etary candidates. Data reduction pipeline, observational strategy, 
and analysis tools were developed to optimize the follow-up of 
transiting candidates. 

Observations were made with the SOPHIE high efficiency 
mode, with a spectral resolution of 39,000 at 550nm, and the 
slow CCD read-out mode. The obj_AB observing mode was 
used, i.e., without acquisition of the simultaneous calibration 
lamp in order to monitor the background sky on the second fiber. 
The typical intrinsic stability of SOPHIE does not require using 
the simultaneous calibration, and the sky fiber is crucial for re- 
moving the scattered moonlight (Santerne et al. 1201 lbi . SOPHIE 
spectra were reduced with the online pipeline. Radial veloci- 
ties were obtained by computing the weighted cross-correlation 
function (CCF) of the spectra using a numerical spectral mask 
corresponding to a G2V star (Baranne et al. 119961 Pepe et al. 
120021 . 

Observations were carried out from 2010 July 15 to 2010 
September 13. A total of 28 measurements with 1-hour expo- 
sure times were done corresponding to an equivalent of about 

1 http://archive.stsci.edu/kepler/data_search/search.php 
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Fig.l. Phase-folded SOPHIE radial velocity of KOI-4T0. The 
three curves correspond to the detection limit at 1, 2, and 3 <x. 



four nights spread over 24 different nights. Such flexibility in the 
schedule of the observations is crucial to properly cover orbital 
phases of transiting candidates. 

4. KOI-410 

We made two SOPHIE measurements on this target. Phase- 
folded radial velocities are shown in Fig.[T]and present no signif- 
icant variation at the level of 120 m s _1 . At 3-cr we can exclude 
a transiting companion with mass greater than 3.4 Mj up . The 
transit light curve, shown in Fig. [2] is however, significantly V- 
shaped. The most probable contaminating star is KIC-5449780, 
which is an 18.4 Kepler magnitude star (36 times fainter) located 
at 9.6 arcsec. We estimated that, if this star was an eclipsing bi- 
nary, it would cause a a displacement of the photocenter motion 
in the X and Y directions of more than 5 milli-pixels during the 
transit, a value excluded by the data (rms of 0.4 milli-pixel in 
both X and Y). This excludes the faint star located at 9.6 arc- 
sec as the origin of the transit. We cannot exclude, however, an 
unresolved diluted eclipsing binary (DEB). To exclude a transit- 
ing giant planet with mass greater than 0.3 Mj up at 3-cr requires 
additional RV measurements at the extrema phases with uncer- 
tainties of 10 m s _I , which is beyond the capabilities of SOPHIE 
for such a faint star. 

5. KOI-552 

We made two SOPHIE measurements on this target (see Fig. [3). 
They show a large radial-velocity variation in phase with the 
Kepler ephemeris. Assuming a circular orbit - which is reason- 
able for a three-day period - the semi-amplitude K is 21.35 + 
0.14 kms~' and reveals a low-mass transiting star with a mass 
ratio of 0.15. 
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Fig. 2. Unbinned phase-folded transit light curve of KOI-410. 
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Fig. 3. Phase-folded SOPHIE radial velocity of KOI-552 



The radius and mass for the central star are estimated to be 
1 .06 R G and 1 . 10 M by Borucki et al. d2011bl ). respectively. The 
derived radius of the companion is then 1.0 Rj up - Such a radius 
seems small for an M-dwarf companion in this range of mass, 
where the typical radius is expected to be 0.17 R Q . Borucki et 
al. (1201 Vo\ recognize that some of the characteristics listed for 
the stars are uncertain, especially surface gravity. This might in- 
troduce errors in stellar diameters up to 25%. Our two SOPHIE 
spectra unfortunately have too low signal-to-noise ratio (S/N) to 
determine the spectral type of KOI-552. From J & K magnitude 
and following the same procedure as for KOI-428b ( Santerne et 
al. 1201 1 at and KOI-423b (see section [6721 in this paper), we esti- 
mate T e g =6560±150K assuming an [Fe/H]=-0.47 as provided in 
the KIC catalog, a value significantly higher than the estimation 
given by Borucki et al. (2011), 7^=60 18K. Assuming a solar 
metallicity changes the estimated r e ff within the uncertainties. If 
we assume that the effective temperature of the star is close to 
6560K, we may expect a larger stellar radius and thus a larger ra- 
dius for the M-dwarf. We also note that the uncertainties on the 
impact parameter are large (b=0.79 + 0.24), because the transit 
is V-shape (see Fig. 0. An underestimation of the impact param- 
eter could also explain the discrepancy in the M-dwarf radius. 
Our two SOPHIE spectra unfortunately have too low S/N to de- 
termine the spectral type of KOI-552. Additional spectroscopy 
measurements are required to better characterize the parameters 
of the detected transiting companion. 

From the MAST database, we downloaded the raw Kepler 
light curve containing photometric data from both the first and 
second quarters. We corrected the Ql and Q2 light curves sep- 
arately for long-term trends due to systematics and removed a 
few of the data points clearly affected by instrumental effects, 
such as those occurring just after the two safe mode events in 
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Fig. 4. A zoom of the KOI-552 light curve (Q2 data) showing 
flux variations due to starspots and, superposed, transits by a 
low-mass stellar companion. 



the Q2 light curve (Haas et al. l20T0l Jenkins et al. [20T0l . We 
estimated the flux excess due to stellar crowding in the Kepler 
mask by comparing the values of the median flux at the be- 
ginning of the raw Ql and Q2 light curves with those of the 
corresponding error-corrected light curves (Jenkins et al. l2010l l. 
Finally, we subtracted the flux excess from each point of the 
Ql and Q2 raw time series, separately. The treated light curve 
presents clear variations due to rotating stellar spots (see Fig.|4j». 
The auto-correlation function of the light curve was computed 
after removing the transits, and a clear peak was found at P, ot = 
3.08+0.15 days indicating a possible spin-orbit synchronisation, 
which is expected for this mass ratio. The v sin /* derived from 
the CCF is estimated to 17.6 + 1.0 kms -1 , which agrees with the 
measured rotational period and a stellar radius greater than 1 R . 
The transit light curve shown in Fig.[5]presents some transits that 
are shallower than others, which may be interpreted by the fact 
that the transiting companion crosses dark stellar spots. The stel- 
lar variability with a period equal or close to the orbital period 
prevents us from properly determining the ellipsoidal modula- 
tion and the beaming effect (van Kerkwijk et al. I20101 >. 

The Kepler Q2 data, released on February 2011 after our 
spectroscopic observations, permit us to detect the secondary 
transit (see Fig . O ■ The secondary transit is clearly too deep to be 
compatible with a perfect reflecting hot Jupiter with the orbital 
parameters of KOI-552 published in Borucki et al. (1201 lbl) . The 
occultation depth of ~ 250 ppm would imply a geometric albedo 
of 2. Assuming a zero albedo, r e ff =6560 K for the primary star, 
and radius ratio R c //?*=0.0967, the effective temperature of the 
companion is estimated to be r e ff=3300 K, which corresponds 
to an M7 dwarf. Although we cannot exclude the possibility that 
a fraction of the observed depth might be due to reflected light, 
the secondary transit indicates a low-mass star companion. 

The V-shape of the transit, the depth of the secondary eclipse, 
and the possible synchronisation between the companion and the 
star provide some evidence that this candidate was a very un- 
likely planetary candidate, although it was ranked as priority 2 
by the Kepler team, . Our two RV measurements definitively es- 
tablish this candidate as a low-mass star. 
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Fig. 5. Phase-folded light curve showing the unbinned primary 
transit (upper panel) and the binned secondary transit (lower 
panel) of KOI-552. The lower pannel also shows the expected 
model of the secondary transit for a hot Jupiter with a geometric 
Albedo of 1 . 



6. KOI-423 

6. 1 . Kepler light curve analysis 

The target KOI-423 was observed by Kepler with a temporal 
sampling of 29.4 min for 122.2 days. The publicly-available 
Kepler light curve contains 5708 photometric measurements in 
total, after 135 data points were discarded by the Kepler pipeline 
as affected by systematic effects (Jenkins et al. 2010). The raw 
Ql and Q2 light curves were downloaded from the MAST and 
treated in the same way as for KOT552. The obtained Kepler 
light curve, shown in Fig. [6j exhibits five transits with a period 
of 21.1 days, a depth slightly shallower than 1% and a duration 
of about six hours. One additional transit should be observable 
but, unfortunately, it fell entirely inside a gap taking place from 
2455014.5 to 2455016.7 (BJD). Furthermore, the fourth transit 
is only partially visible because of another data gap just before 
the ingress. 
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Fig. 6. The Kepler light curve of KOI423 with a temporal sam- 
pling of 29.4 min showing five transits and variations due to stel- 
lar activity. 




Fig. 7. Phase-folded radial velocity of KOI-423 with the best 
Keplerian orbit fit. Open dots correspond to measurements af- 
fected and corrected by the scattered moonlight. 



The KOI-423 light curve shows flux variations due to the 
presence of active regions on the photosphere, whose visibility 
is modulated by the stellar rotation. The peak-to-peak amplitude 
is ~ 0.6%. 

The standard deviation of the light curve, computed in a ro- 
bust way after removing low frequency variations with a sliding 
median filter, is 2.5 10~ 4 , which is compatible with the median of 
the errors of the single photometric measurements, i.e. 2.1 10~ 4 . 

To estimate the rotation period of the star, we computed the 
auto-correlation function of the Kepler light curve, after remov- 
ing the five transits, and found a clear peak at f rot = 4.35 + 
0.21 days. We also computed the Lomb-Scargle periodogram of 
the light curve (Scargle et al. 1 19821 1. and found power at 4.35, 
2.17, and 1.45 days corresponding to the stellar rotation period 
and its two first harmonics. 



6.2. SOPHIE spectra and radial velocities 

Thirteen measurements were made on this candidate with 
SOPHIE from 2010 July 26 to 2010 September 10. They are 
listed in Table [2] Seven of them were strongly affected by the 
scattered moonlight but corrected for radial velocity analysis us- 
ing the sky spectra acquired simultaneously. In some exposures, 
the S/N on the sky fiber B, which points on the sky background, 
was even comparable to or greater than on the target fiber A, 
which is located on the target. We conservatively doubled the 
uncertainties of these seven measurements and did not use them 
to compute the bisector span. 

The radial velocities shown in Fig. [7] present a clear varia- 
tion in phase with the Kepler ephemeris and with a slight but 
significant eccentric orbit. The bisector span, listed in Table [2] 
and shown in Fig. [8] does not reveal any significant variations 
at a level more than ten times smaller than the radial velocity 
variations allowing all scenarios of diluted and blended eclips- 
ing binaries to be excluded. 

We performed the spectroscopic analysis of the parent star 
using SOPHIE spectra. Individual spectra are too poor in qual- 
ity to allow a proper spectral analysis. The six SOPHIE spec- 
tra not affected by the scattered moonlight were shifted to the 
barycentric rest frame and co-added order per order. We finally 
got a co-added spectrum with an S/N of about 110 per pixel on 
the continuum at 550 nm. From the analysis of a set of isolated 
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lines, we derived a vsinz* of 16+2.5 kms -1 . The spectroscopic 
analysis was carried out using the same methodology as for the 
CoRoT planets described in detail in Bruntt et al. ( 120101 1. 

We limited the abundance analysis to the iron ions, because 
of the high vsim* of the star and the quite moderate spectral 
resolution and S/N of the combined spectrum. In total 53 Fel 
and 11 Fell spectral lines were used for the temperature and 
logg determination. Following the Bruntt et al. (2008) method- 
ology, abundances were measured differentially with respect to a 
SOPHIE solar spectrum. We found a marked underabundance of 
iron, with [Fe/H]=-0.29 ±0.1. This result was further checked 
using the SOPHIE CCFs. Following the relations established by 
Boisse et al. (2010), we derived a vsin/* of 12.6 + 1 kms -1 
and an [Fe/H] of -0.18 + 0.1 dex. The metallicity is l-cr higher 
than the spectroscopic determination. The vsin/* value esti- 
mated from the CCF is 3-cr lower than then spectroscopic de- 
termination, which may come from the fact that the relation was 
calibrated assuming solar metallicity. The rotation velocity of 
KOI-423 derived from the stellar radius and the rotational pe- 
riod estimated from the light curve is 16.2 +1.7 kms 1 , in full 
agreement with the spectroscopic determination, indicating that 
i+ may be close to 90 deg. 

The temperature derived from the Fe i - Fe n equilibrium is 
T eff = 6260 + 140 K. This is 2 cr higher than the value of 5992 K 
quoted for the star in the Kepler Input Catalog. From infrared 
magnitudes J and K given in the 2MASS catalog, a redden- 
ing of E(J-K)=0.0 91 (fro m Cardelli et al. [1989] ), and Eq. 3 in 
Casagrande et al. (2010) assuming [Fe/H] = -0.3 dex, we de- 
rived an effective temperature T e ff = 6380 ± 150 K. This esti- 
mate agrees with our spectroscopic determination. As already 
mentioned for KOI-428 (Santerne et al. 1201 lab , the infrared de- 
termination of the effective temperature from Casagrande et al. 
d2010l l. based on stars of luminosity class IV and V, seems to be 
more reliable than the one used for the Kepler Input Catalog. 

The surface gravity was estimated using pressure-sensitive 
lines: the Mgilb lines, the Nai D doublet, and the Cai at 6122 
and 6162 A. As already noted by Bruntt et al. (120101) . the results 
obtained on the Mg i triplet vary from one line to another due 
to the difficulty in the normalization of these broad lines. Fitting 
each of the observed lines of the triplet with synthetic spectra 
calculated using MARCS models, we found log g = 3.90 + 0.25. 
The Na i D lines gave log g = 4.0 ± 0.2 and the Ca i lines log g = 
4.1 ± 0.2. The latter are in good agreement with the value of log 
g = 4.2 obtained with VWA from the agreement between the Fe i 
and Fe n abundances. Considering the issue on the normalization 
for broad lines such as Mgi, we adopted as a final value the 
gravity derived from the Ca i pressure sensitive lines, that is, log 
g = 4.1 ±0.20. 

6.3. System parameters 

System parameters were derived by performing the transit mod- 
eling and the Keplerian fit of the radial velocity measurements 
simultaneously. Transit fitting was carried out following the for- 
malism of Gimenez (2006 2009) after removing stellar variabil- 
ity in the out-of-transit light curve to correctly normalize the 
transits. For this purpose, we fitted a third-degree polynomial 
to the 9h intervals of the light curve before the ingress and after 
the egress of each transit. We discarded the fourth transit as only 
half visible, which prevented us from fitting the polynomial be- 
fore the transit ingress and thus normalizing it as for the other 
ones. 



Table 2. SOPHIE measurements of KOI-423 



BJD 


RV 




BIS 


S/N/pix 


-2455 000 


[kms- 1 ] 


[kms" 1 ] 


[kms" 1 ] 


@550nm 


403.49308" 


-0.173 


0.225 


- 


8.2 


425.43895 


-0.059 


0.056 


-0.040 


18.2 


429.45136" 


0.877 


0.092 


- 


21.1 


431.39644" 


1.104 


0.104 


- 


20.3 


432.33993" 


1.440 


0.150 


- 


16.4 


433.47474" 


0.952 


0.127 




16.5 


435.45034" 


0.598 


0.168 




11.9 


438.40205" 


-0.765 


0.167 




11.0 


441.36981 


-1.324 


0.041 


-0.058 


20.7 


442.39937 


-1.278 


0.056 


-0.088 


18.4 


443.37619 


-1.145 


0.037 


-0.009 


21.9 


445.38001 


-0.528 


0.043 


0-0.198 


22.7 


450.39197 


0.754 


0.037 


0-0.145 


23.5 



" Measurements corrected from the scattered moonlight. 



KOI-423 SOPHIE 




RV |km.B 'J 



Fig. 8. Bisector span versus radial velocity of KOI-423 showing 
no correlation. 



The nine free parameters of our global best fit are the orbital 
period P, the transit epoch T u , the transit duration d tr , the ratio 
of the companion to stellar radii R c /R+, the inclination i between 
the orbital plane and the plane of the sky, the Lagrangian orbital 
elements h = e sin to and k = e cos u>, where e is the eccentricity 
and to the argument of the periastron, the radial-velocity semi- 
amplitude K, and the systemic velocity Vo. The two nonlinear 
limb-darkening coefficients u+ — u a + u\, and m_ = u a - ut,^ 
were fixed in the transit modelling. They cannot be derived as 
free parameters because the transit ingress and egress are not 
sampled well, given the coarse temporal sampling of the Kepler 
light curves. The adopted limb-darkening coefficients u a and uu 
for the Kepler bandpass were taken from Sing's 120101 tablefl 
after linearly interpolating at the T e ff, logg, and metallicity of 
the star: u a = 0.303 ± 0.014 and u b = 0.308 ± 0.005, which give 
u+ = 0.611 ± 0.015 and U- = -0.004 ± 0.015. 

The best-fit parameters were found by using the algorithm 
AMOEBA (Press et al. 1992) and changing the initial values of 
the parameters with a Monte-Carlo method to find the global 
minimum of the x 1 - Following Kipping & Bakos d201 11 1 and 

2 u a and ut are the coefficients of the limb-darkening quadratic law: 
IQx)/I(l) = 1 —u a (\ -/i) — u b (\ -fi) 2 , where 1(1) is the specific intensity 
at the center of the disk and /j = cos y, with y the angle between the 
surface normal and the line of sight 

3 http://vega.lpl.arizona.edu/singd/David_Sing/Limb_Darkening.html 
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Santerne et al. (1201 lal l. the transit modeling was performed with 
a temporal sampling five times denser than that of Kepler, i.e. 
one point every 5.88 min, and then binning the model samples 
to match the Kepler sampling rate. In such a way, the solu- 
tion of the transit parameters is more accurate than a simple y 2 
minimization between the Kepler measurements and the tran- 
sit model. Indeed, the Kepler sampling of 29.4 min inevitably 
smears all sharp changes occurring, such as transit ingress and 
egress (Kipping & Bakos 201J}. 



Table 3. KOI-423 system parameters. 
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Fig. 9. Top panels: Unbinned phase-folded Kepler transit light 
curve of KOI-423. The red solid line shows our oversampled 
transit model. Bottom panels: the residuals from the best-fit 
model. 

Fitted and derived system parameters are listed in Table [3] 
together with their 1-cr errors estimated using a bootstrap pro- 
cedure. The last works as follows: 1) subtracts the best solution 
to the data; 2) shifts the photometric residuals, thus taking also 
possible correlated noise into account; 3) shuffles a fraction of 
the radial velocity residuals (typically 1/e ^ 37%); 4) adds the 
subtracted solution; and 5) again carries out the global best fit. 
The limb-darkening parameters u+ and U- in the transit mod- 
eling were allowed to vary within their error bars related to the 
atmospheric parameter uncertainties. More than one thousand it- 
erations were used for this procedure. FigureQshows the phase- 
folded radial velocity of KOI-423 with the best Keplerian fit. 
Figure [9] shows the phase-folded transit light curve of KOI-423 
and, superposed, the transit model with a sampling rate of 5.88 
min. 

The star's mass and radius were estimated using a grid of 
STAREVOL evolutionary stellar tracks calculated for the CoRoT 
exoplanet program (Palacios, 2007, private com). We used the 
regular approach for exoplanet host stars, which is the compar- 
ison of the star location in a (T e fr, M l J 3 /R+) HR-plane to stel- 
lar evolutionary tracks with the appropriate range of metallicity. 
Using the M^ 3 //?*. from the transit modeling (see Table [3) and 
the effective temperature, we estimated the mass M* = 1.10 
from the tracks and deduced the radius from this mass value and 
the / Revalues. Combining these mass and radius values we 
get log g = 4.20 ± 0.2, which is in good agreement with the spec- 
troscopic value. With the adopted stellar parameters, we derived 
for the transiting companion M c = 18.0 + 0.92 Mj up and R c = 
1.22 ± 0.11 R Jup . 

The distance of the star can be estimated from the V magni- 
tude (mv = 14.46), a bolometric correction BC= -0.03 for F7IV 
stars (Straizys & Kuriliene 1981), the adapted solar bolometric 



Parameters 



Value 



Transit epoch T„. [BJD] 
Orbital period P [days] 
Transit duration T 14 [h] 
Radius ratio R c /R+ 
Inclination i [deg] 
h = e sin u> 
k = e cos a* 

Semi-amplitude K [kms~'] 
Systemic velocity V u [kms~'] 

Orbital eccentricity e 
Periastron argument o> [deg] 
Scaled semi-major axis a/R* 
Ml /3 /R* [solar units] 
Stellar density p+ [g cnr 3 ] 
Impact parameter b 

Effective temperature T e ff [K] 
Surface gravity log g [dex] 
Metallicity [Fe/H] [dex] 
Rotational velocity v sin [ km s _I ] 
Spectral type 

Star mass [M ] 
Star radius [R Q ] 
Distance of the system [pc] 

Orbital semi-major axis a [AU] 
Companion mass M c [Mj up ] 
Companion radius ^ c [Rj up ] 
Companion density p c [g arT 3 ] 
Equilibrium temperature T cq [K] " 



Jt7li.J7J5 00|)05 

21.0874 + 0.0002 
6.02 ± 0.09 
0896 +000U 
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23.8 + J : = 
740 +0 057 
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-0.18 

6260±140 
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-0.29±0.10 

16+2.5 

F7IV 

j 1 +o.o7 

1,J -o.io 
1200±250 

0.155 ±0.003 

18.00 + "'-? 



i 99+0.12 
I - -0.10 

12.40 + ^^ 
905 + ^ 



" black body equilibrium temperature for an isotropic 
planetary emission. 



magnitude (Mboi,o=4.74) recommended by Torres (2010J, the 
reddening coefficient Ay =0.423 provided by MAST database, 
and the stellar luminosity derived from the stellar radius and T e ff. 
We find d- 1200+250 pc. The age was estimated to be 5.1 + 1.5 
Gyr and 4.4+1.0 Gyr from the STAREVOL and the CESAM 
stellar evolution models, respectively. All the derived stellar pa- 
rameters are reported in Table [3] 

6.4. KOI-423b and the evolution of brown dwarfs and 
massive exoplanets 

Figure [TO] compares the size of KOI-423b to other known tran- 
siting objects from 5 Mj up to the M-dwarf mass range. Its size 
is large, and it appears to be inflated, although comparable to 
other massive planets such as XO-3b (Johns-Krull et al. 2008). 
However, given its high mass, it lies farther from "classical" evo- 
lution tracks, such as those for isolated objects from Baraffe et 
al. d2003), than other known objects of this class. The problem 
is made even more acute when one considers that KOI-423b is 
relatively cold (r eq = 905 K), in contrast to most other objects in 
this plot, implying that irradiation effects cannot be held respon- 
sible for its anomalously large size. 

In Fig. QT| we compare the constraints in age and radius 
obtained from photometric, spectroscopic, and CESAM stel- 
lar evolution models to dedicated evolution models including 
irradiation (see Guillot & Havel I20TT1 and Havel et al. I20TT1 
for a description of the method). We first calculate a "stan- 
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Fig. 10. Mass-radius diagram of massive transiting planets, 
eclipsing brown-dwarfs, and low-mass stars. Open symbols cor- 
respond to host stars with T e ff < 6000 K. Dotted lines correspond 
to theoretical isochrones with ages 0.5, 1, and 5 Gyr (from top 
to bottom) from Baraffe et al. (20Q3). 
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Fig. 11. Derived constraints on the radius of KOI-423b as a 
function of age compared to evolution models. The dots rep- 
resent solutions within lcr (red), 2<x (orange), and 3cr of the 
inferred (r e ff,p*) obtained from stellar evolution models. From 
bottom to top, the lines correspond to evolution calculations for a 
hydrogen-helium brown dwarf with Y = 0.25 using a "standard" 
evolution model (black line), a model in which the greenhouse 
factor has been arbitrarily multiplied by 10 (purple line), and a 
model in which the atmospheric opacities have been arbitrarily 
multiplied by 10 (blue line). 



dard" evolution track with a simplified semi-gray atmospheric 
boundary condition that has been parametrized to agree with 
radiative transfer calculations and evolution calculations (see 
Guillot 120 10) . Specifically, we adopt a thermal (infrared) opac- 
ity coefficient k± = 0.04 cm 2 g~' and a visible opacity k v - 
0.024 cm 2 g -1 . The corresponding evolution track is found to be 
at the 3<x lower limit of the main-sequence solutions. When we 
increase the greenhouse effect by a factor 10 by decreasing the 
visible opacity by the same factor, the consequence on the evolu- 
tion is small. This is because the deep atmosphere is convective 
and weakly affected by a deeper penetration of the stellar light. 
Similarly, recipes commonly adopted to explain the large size 



of exoplanets of smaller sizes, including tides, dissipation of ki- 
netic energy, or increased interior opacities (see Guillot 2008) 
have little effect on this very large planet. Even the possibility 
that only 50% of the atmosphere would be allowed to transfer 
heat efficiently (proposed by Bouchy et al. 1201 li as one possibil- 
ity to explain the large size of CoRoT-15b) does not work in this 
case. We also checked that the energy dissipation from the circu- 
larization of the orbit by tides is quite inefficient at the distance 
of KOT423b.To obtain the required dissipation rates, one would 
need a tidal quality factor Q p of about 10, i.e., at least 3 orders of 
magnitude over Jupiter's. In fact, the only way we found to rec- 
oncile models and observations within ~ lcr is to invoke an ad- 
hoc increase in the atmospheric opacities, by an order of magni- 
tude or more (we used /c t h = 0.4 cm 2 g -1 and k v = 0.24 cm 2 g~ l ). 
This is in essence similar to what is proposed by Burrows et al. 
(2007) for planets with lower masses. This yields a slower cool- 
ing of the object that remains large for a longer time. 

However, the comparison between KOI-423b and CoRoT- 
3b (Deleuil et al. 2008) adds to the puzzle: the two compan- 
ions have different sizes (CoRoT-3b lies where theory predicts 
it should be, KOT423b is significantly larger), but they other- 
wise have very similar characteristics (mass, characteristics of 
the parent star). KOT423 has vsin/ = 16 + 2.5 km s -1 , whereas 
CoRoT-3 has v sin i = 17 ± 1 km s _1 , showing that both are rel- 
atively fast rotators. CoRoT-3 has a metallicity close to the Sun 
([M/H] = -0.02 ± 0.06) when KOI-423b is metal-poor. The 
smallest massive planet HAT-P-20b (Bakos et al. 1201 Oi l orbits a 
host star with the strongest metallicity index [Fe/H] = 0.35. On 
one hand, this follows the trend observed at lower masses that 
planets around metal-rich stars contain more heavy elements and 
are therefore on average smaller (see Guillot 2008, Bouchy et al. 
2010). On the other hand, one would expect that the atmosphere 
of KOI-423b is poorer in heavy elements, hence yielding smaller 
opacities, and potentially, a smaller size. 

One natural possibility would be to invoke age as a way to 
distinguish between the two objects. However, given the esti- 
mated relatively high age, it is difficult to imagine that KOI- 
423 would be on the pre-main sequence (see Fig. fTTb . whereas 
CoRoT-3b would be much older. We cannot pretend that we un- 
derstand the large radius of KOI-423. 

6.5. Comparison to the population of massive exoplanets 

For these systems with a known substellar companion, we plot in 
Fig.Q~2]the stellar effective temperature as a function of the com- 
panion mass for transiting objects or M sin i value for nontran- 
siting objects. The data for low-mass companions (< 13Mj up ) 
was obtained from exoplanets.org (see Wright et al. 2011} and 
exoplanet.eu. For objects with higher masses, no publicly avail- 
able database exists, and we manually extracted parameters for 
the following transiting substellar systems: OGLE-TR-106 (Pont 
et al. |2005al . O GLE-T R-122 (Pont et al. I2005bl l. OG LE-TR - 
123 (Pont et al. 12006) . H AT-TR- 205-013 (Beatty et al. 12007) . 
CoRo T-15 (Bouchy et al. 12017) . W ASP-30 (Anderson et al. 
12017) , and LHS6343C (Johnson et al. 12077) . And we did it for 
non-transiting short-period (P < 30 days) substellar systems: 
NLTT 41135C (Irwin et al. l2070) . H D1893 10 (Sahlmann et al. 
12077) . and HD180777 (Galland et al. 12006) . 

Although our ensemble of objects is clearly incomplete and 
derived from inhomogeneous databases, the lack of short period 
(P < 10 days), massive (M c > 5Mj up ) companions around 
G-type stars (T e ff between 5000 and 6000K) appears significant 
and puzzling. In comparison, these massive companions are rou- 
tinely found next to F-type stars (including KOI-423). Given that 
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Fig. 12. Values of the stellar effective temperature and mass of 
the companion for systems containing a companion in the planet, 
brown dwarf, or M-dwarf mass range. Transiting/eclipsing sys- 
tems are shown by circles, and systems only detected in radial 
velocimetry are shown by trianges (in which case the mass value 
corresponds to the mass times the sine of the inclination of the 
orbit with respect to the observer). The size (area) of the sym- 
bols is inversely proportional to the orbital period (see labels). 
The metallicity of the parent star is shown by the symbol color 
(see labels). Some objects deemed important for understanding 
the link between these parameters are labeled. The error bars 
on the masses are small (~ 10%) and would be barely visible 
in this diagram. The error bars on the effective temperatures are 
believed to range from ~ 100K to up to ~ 300 K (for example 
for the faint OGLE-TR-122). 



selection biases would tend to favor their detection around G 
dwarfs, we believe that their absence is real. 

The empty region is bounded very approximatively by the 
following objects: HAT-P-20b, a massive 7.5 Mj up , anomalously 
small (see Fig. ITOb planet orbiting a K7 dwarf in 2.88 days; 
CoRoT-2b, an anomalously large 3.3 Mj up planet orbiting a G7 
dwarf in 1.74 days; CoRoT-14b, a 7.6Mj llp planet orbiting an 
F9V star in 1.5 days; KOI-423b; WASP-30b, a 60M Jup brown- 
dwarf orbiting a relatively young (1-2 Ga) F8V star. At the top 
of the diagram, OGLE-TR-122b is a 96 M Jup M dwarf on a 7.27- 
day orbit around an apparent G dwarf. However, in this case, 
given the faintness of the star and its poorly-determined r e ff 
(5700 + 300 K), it may in fact also belong to the F-dwarf cat- 
egory. 

Although this is not shown in Fig. Q2] we believe that this 
trend continues well into the M-dwarf regime. For example, 
a study of 18 OGLE transit candidates with periods between 
0.8 and 13.9 days (Bouchy et al. 2005) identified two planets 
(OGLE-TR-lOb and OGLE-TR-56b; present in Fig. [12]), 7 bina- 
ries consisting of an M-dwarf orbiting an F-dwarf and 1 M-dwarf 
around G-dwarf binary. 

It is unlikely that massive companions do not form close to 
G-dwarfs or do not migrate there: the difference in mass be- 
tween G-dwarfs and F-dwarfs ( 20%) is small compared to the 
extension in mass of the "G-dwarf companions desert" (from 
about 10Mj llp to well within the stellar-mass regime). Following 
Bouchy et al (1201 11 1, we propose that close-in massive planets, 
brown dwarfs, or M-dwarfs may survive when orbiting orbiting 
close to F-dwarfs, but not when close to G-dwarfs. This could be 



because around G-dwarfs, the star would rapidly spin down and 
tidal interactions would lead to a migration and eventual engulf- 
ment of its companion. In contrast, around F-dwarfs, because of 
their small or absent convective zone, the much weaker braking 
(see e.g. Barker & Ogilvie 2009) will avoid the loss of angular 
momentum and the rapid decay of the companion's orbit. 

Further work is required including both a proper statistical 
analysis of the observed correlations and models that combine 
the formation context, tidal, and magnetic interactions between 
the companions. Kepler will certainly help in providing a much 
more well-sampled comparison between F and G-dwarfs giant 
planetary companions. 



7. Conclusions 

The SOPHIE follow-up of four Kepler transiting candidates, se- 
lected with the aim of privileging massive planets and brown 
dwarfs, revealed at least three real massive transiting compan- 
ions with Jupiter sizes: KOI-428b, a 2.2 Mj up transiting planet 
(Santerne et al. 1201 lab : KOI-552b, a ~ 0.17 M Q eclipsing low- 
mass star; and KOI-423b, a 18.0 Mj up transiting companion. The 
status of the fourth candidate KOI-410 has not yet been solved, 
but we can exclude a companion heavier than 3.4 Mj up , and we 
still may suspect a possible background eclipsing binary. Our 
result shows that at least three Kepler candidates over the four 
selected are real Jupiter-size transiting objects, but only one of 
them, KOI-428b is in the mass domain of Jupiter-like planet. 

KOI-423b is one of the rare types of transiting companions, 
like CoRoT-3b (Deleuil et al.. 120081 ) that is located in the over- 
lapping region between massive planet and brown dwarf and 
which may be considered either as a "superplanet" or a low- 
mass brown dwarf. As discussed in section l6~4l its anomalous 
bloated radius is clearly not understood. Additional high S/N 
and high-resolution spectroscopic observations, as well as ad- 
ditional Kepler data, will permit us to better constrain the size of 
KOI-423b. 

From RV surveys, only five massive planets or brown dwarfs 
(M c .sin/ > 5 Mj up ) with short periods (P < 30 days) were 
found including HD162020b (14.4 M Jup ), HD41004Bb (18.4 
M Jup ), HD180777b (25 M Jup ), HD189310b (25.6 M Jup ), and 
NLTT41135C (33.7 M Jup ) respectively orbiting K2V, M2V, 
A9V, K2V, and M5V stars. Except for HD 180777b fou nd in a 
dedicated RV surveys around A-F type star (Galland et al. 2006), 
this lack of short-period massive planets and brown dwarfs prob- 
ably comes from the bias in the selection of RV surveys focusing 
on spectral types later than F5 and on slow rotators. 

After our first SOPHIE campaign, Borucki et al. d201 lbb 
published an update of the Kepler exoplanet candidates totalling 
1235 candidates. Using the same criteria of selection as pre- 
sented in section|2](R > 0.8 R Jup , mv < 14.7 and T eff > 5900 K), 
we found 17 additional candidate^. Radial velocity follow-up 
of these candidates may help to increase the statistic, and to fur- 
ther reveal the properties of massive planets and brown dwarfs, 
as well as to assess the true fraction of false positives in order to 
correctly interpret the Kepler candidate population. 

Acknowledgements. We thank the technical team at the Observatoire de Haute- 
Provence for their support with the SOPHIE instrument and the 1 .93-m telescope 
and, in particular, for the essential work of the night assistants. Financial support 
from the Programme national de planetologie (PNP) of the CNRS/INSU, France, 



4 This list includes the recently announced Kepler- 14b (Bucchave et 
al. 1201 It . a close-in massive planets of 8.4 Mj up orbiting close to an 
F-type star. 



F. Bouchy et al.: KOI-423b : an 18 Mj up transiting companion around an F7IV star 



9 



and from the Swiss National Science Foundation (FNSRS) are gratefully ac- 
knowledged. We also acknowledge support from the French National Research 
Agency (ANR-08-JCJC-0102-01). ASB is supported by a CNES grant. AE is 
supported by a fellowship for advanced researchers from the Swiss National 
Science Foundation. 



References 

Anderson, D.R., Collier Cameron, A., Hellier, C, et al, 2011, ApJ, 726, L19 
Bakos, G.A., Hartman, J., Torres, G., et al., 2010, ApJ, submitted 

larXiv:1008.3388l 
Baraffe, I., Chabrier, G., Barman, T.S., et al., 2003, A&A, 402, 701 
Baranne, A., Queloz, D., Mayor, M., et al., 1996, A&AS, 119, 373 
Barker ,A., & Ogilvie, G.I., 2009, in Cosmic Magnetic Fields: From Planets, to 

Stars and Galaxies, Proceddings of the International Astronomical Union, 

IAU Symposium, 259, 295 
Batalha, N.M., Borucki, W.J., Bryson, S.T., et al., 2011, ApJ, 729, 27 
Beatty, T.G, Fernandez, J.M., Latham, D.W., et al., 2007, ApJ, 663, 573 
Boisse, I., Eggenberger, A., Santos, N.C., et al., 2010, A&A, 523, 88 
Borucki, W.J., Koch, D.G., Basri, G, et al., 2011a, ApJ, 728, 117 
Borucki, W.J., Koch, D.G, Basri, G, et al., 2011b, ApJ, 736, 19 
Bouchy, F, Pont, F, Melo, C, et al., 2005, A&A, 431, 1105 
Bouchy, F, Hebrard, G, Udry, S., et al., 2009, A&A, 505, 853 
Bouchy, F, Hebb, L., Skillen, I., 2010, A&A, 519, 98 
Bouchy, F, Deleuil, M., Guillot, T., et al., 201 1, A&A, 525, 68 
Bruntt, H., De Cat, P., & Aerts, C, 2008, A&A, 478, 487 
Bruntt, H., Deleuil, M., Fridlund, M., et al., 2010, A&A, 519, 51 
Buccha ve, L.A., Latham, D.W., Carter, J.A., et al., 2011, ApJ, submitted 

(arXiv: 1106.55 lOt 
Burrows, A., Hubeny, I., Budaj, J., Hubbard, W. B., 2007, ApJ, 661, 502 
Cardelli, J.A., Calyton, G.C., & Mathis, J.S., 1989, ApJ, 345, 245 
Casagrande, L., Ramirez, I., Melendez, J., et al, 2010, A&A, 512, 54 
Deleuil, M., Deeg, H., Alonso, R., et al., 2008, A&A, 491, 889 
Ehrenreich, D., Lagrange, A.M., Bouchy, F, et al, 2011, A&A, 525, 85 
Galland, F, Lagrange, A.M., Udry, S., et al, 2006, A&A, 452, 709 
Gimenez, A. 2006, A&A, 450, 1231 

Gimenez, A. 2009, in The Eighth Pacific Rim Conference on Stellar 

Astrophysics: A Tribute to Kam Ching Leung, Eds. B. Soonthornthum, S. 

Komonjinda, K.S. Cheng, and K.C. Leung, ASP Conf Series, Vol. 404, 291 
Guillot, T, 2008, PhST, 130, 4023 
Guillot, T., 2010, A&A, 520, 27 
Guillot, T., & Havel, M., 2011, A&A, 527, 20 
Haas, M., Batalha, N., Bryson, S., et al., 2010, ApJ, 713, 1 15 
Havel, M., Guillot, T. Valencia, D., & Crida, A., 2011, A&A, 531, 3 
Irwin, J., Buchhave, L., Berta, Z.K, et al., 2010, ApJ, 718, 1353 
Jenkins, J.M., Caldwell, DA., Chandrasekaran, H, et al., 2010, ApJ, 713, L120 
Johnson, J.A., Apps, K, Gazak, J.Z., et al., 201 1, ApJ, 730, 79 
Johns-Krull, C, McCullogh, P., Burke, C, et al., 2008, ApJ , 677 , 657 
Kipping, D.M., & Bakos, G.A., 2011, ApJ, 733, 36 
Koch, D.G, Borucki, W.J., Rowe, J.F., et al., 2010, ApJ, 713, 131 
Nordstrom, B., Stefanik, R.P., Latham, D.W., et al., 1997, A&ASS, 126, 21 
Pepe, E, Mayor, M. Galland, F, et al, 2002, A&A, 388, 632 
Perruchot, S., Kohler, D., Bouchy, F., et al., 2008, in Ground-based and Airborn 

Instrumentation for Astronomy II, Edited by McLean, I.S., Casali, M.M., 

Proceedings of the SPIE, vol. 7014, 70140J 
Pont, F, Bouchy, E, Melo, C.M.F., et al., 2005, A&A, 438, 1 123 
Pont, E, Melo, C.H.F., Bouchy, E, et al., 2005, A&A, 433, L21 
Pont, E, Moutou, Bouchy, E, et al, 2005, A&A, 433, L21 
Press, W. H, Teukolsky, S. A., Vetterling, W. T, Flannery, B. P. 1992, Numerical 

recipes in FORTRAN, The art of scientific computing (Cambridge: 

University Press), 2nd edn 
Sahlmann, J., Segransan, D., Queloz, D., et al., 2011, A&A, 525, 95 
Santerne, A., Diaz, R E, Bouchy, E, et al, 2011, A&A, 528, 63 
Santerne, A., Endl, M., Hatzes, A., Bouchy, E, Moutou, C, et al., 

2011, in Detection and Dynamics of Transiting Exoplanets, St. Michel 

l'Observatoire, France, Eds F. Bouchy; R. Diaz; C. Moutou; EPJ Web of 

Conferences, vol 11, 02001 
Scargle, J. D. 1982, ApJ, 263, 835 
Sing, D. K. 2010, A&A, 510, A21 
Straizys, V, & Kuriliene, G, 1981, Ap&SS, 80, 353 
Torres, G, 2010, AJ, 140, 1158 

Udry, S., Mayor, M., Naef, D., et al, 2002, A&A, 390, 267 

van Kerkwijk, M.H., Rappaport, S.A., Breton, R.P., et al., 2010, ApJ, 715, 51 

Wright, J.T., Fakhouri, O., Marcy, G.W., et al., 2011, PASP, 123, 412 



